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Abstract 

Within Density Functional Theory, we have calculated the energy of the transitions from the 
ground state to the first two excited states in the electron bubbles in liquid helium at pressures from 
zero to about the solidification pressure. For 4 He at low temperatures, our results are in very good 
agreement with infrared absorption experiments. Above a temperature of ~ 2 K, we overestimate 
the energy of the Is — lp transition. We attribute this to the break down of the Franck-Condon 
principle due to the presence of helium vapor inside the bubble. Our results indicate that the 
Is — 2p transition energies are sensitive not only to the size of the electron bubble, but also to its 
surface thickness. We also present results for the infrared transitions in the case of liquid 3 He, for 
which we lack of experimental data. 

PACS numbers: 67.40.Yv , 67.55.Lf , 68.03.Cd , 47.55.Dz 
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I. INTRODUCTION 



Excess electrons in liquid helium are known to form electron bubbles. This is so because 
of the strong repulsion between the helium atoms, which are very weakly polarizable, and the 
intruder electrons. Experimental and theoretical spectroscopic studies on electron bubbles 
have been carried out for many years since the pioneering works of Northby and Sanders, 1 
and of Dexter and coworkers.—^ Only some twenty years later, Grimes and Adams were 
able to carried out a detailed analysis of the infrared-absorption spectrum of the electron 
bubble in liquid 4 He. In a first work, 4 they have used a photoconductive mechanism to 
detect the transitions which operated only over a limited region of the pressure-temperature 
(P — T) plane because it appeared to be associated with trapping of bubbles on vortices 
in the superfluid. Later on, they have improved their experimental apparatus and have 
observed the electronic transitions in direct infrared absorption. 5 This has allowed them to 
study the spectrum of electron bubbles in a wider region of the P — T plane from P = 
to the solidification pressure, and up to temperatures above 4 K instead of the maximum 
T ~ 1.6 K reached in photocurrent experiments. 

Recent experiments on cavitation in liquid helium have renewed the interest in the study 
of single electron bubbles.^£ Multielectron bubbles have also prompted some theoretical 
activity j*^ - while the stability of multielectron bubbles is not yet clarified , 10 ! 11 single electron 
bubbles are stable entities for pressures above saturation pressure. Electron bubbles in 
helium droplets are metastable objects whose lifetime has been measured for both helium 
isotopes. 12 

The simplest model to address electron bubbles in liquid helium supposes that they are 
confined in a spherical well potential of radius R. The total energy of the electron-helium 
(e-He) system is then written clS cL function of R 

U(R) = E e + 4nR 2 (r + -ttR 3 P , (1) 

3 

where E e is the ground state electronic energy, P is the pressure applied to the system, 
and a is the surface tension of the liquid. For an infinite spherical well potential, E e = 
ir 2 h 2 1 (2m e R 2 ) . This model is able to qualitatively reproduce the experimental infrared- 
absorption energies, and is simple enough to allow to calculate shape fluctuations of electron 
bubbles and line shapes.— It can be refine d 2 ^ 4 ' 5 by taking into account the finite depth of 



the well Vq, which is about 1 eV, Ref. 
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Once U(R) has been minimized with respect to 



R and the radius of the equilibrium bubble R eq has been determined, it is easy to obtain 
the energies of the (n, I) excited states and to compute the transition energies to the ground 
state if the Franck-Condon principle holds, i.e., if the absorption of a photon by the electron 
bubble is a process much faster than the time needed for the helium bubble to adapt itself 
to the electron wave function in the excited state. 

Another key ingredient entering Eq. (JIJ is the surface tension a, which for 4 He at zero T 
and P is about 0.274 K A~ 2 , Ref. [l5j This yields R eq = 18.9 A. The surface tension is only 
known along the liquid-vapor coexistence line. Consequently, to use Eq. one has to rely 
on model calculations of a(P,T). The situation has been discussed in detail by Grimes and 
Adams, who concluded that in order to perfectly fit their experimental results with Eq. (Q), 
one has to take the depth of the well Vq P-dependent and a surface tension a independent 
of P, which seems difficult to justify. For instance, early calculations*^ indicated that the 
surface tension a nearly doubles from P = to 25 atm. For the sake of completeness, we 
present in the Appendix a calculation of <j(P) that also yields this increase, although we will 
not make a substantial use of these results in our study. It is worthwhile to stress that when 
both Vo and a were made P-dependent following the expected pressure dependences, the 
agreement between theory^ and experiment was only qualitative, see Fig. 2 of Ref. |4j. For 
this reason, these authors concluded that more sophisticated calculations of the properties 
of electron bubbles were called for, incorporating the fact that the bubble density profile is 
not abrupt but have a finite surface thickness. 

In this work we present a theoretical description of infrared-absorption Is — lp and 
Is — 2p excitation energies of electron bubbles in liquid helium. It belongs to the class 
of the more sophisticated calculations suggested by Grimes and Adams, since it is not based 
on ad-hoc assumptions, but on the use of a finite-temperature density functional (DF) 
approach in conjunction with a realistic electron-helium effective potential. DF methods 
have become increasingly popular in recent years as useful computational tools to study 
the properties of classical and quantum inhomogeneous fluids, 18 especially for large systems 
for which these methods provide a good compromise between accuracy and computational 
cost, yielding results in agreement with experiment or with more microscopic approaches. 
In the frame of DF theory, the properties of an electron bubble approaching the surface of 
liquid 4 He have been studied by Ancilotto and Toigc 19 using the so-called Orsay-Paris (OP) 
zero temperature finite-range DF^ and the pseudopotential proposed in Ref. |2l| as e-He 
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interaction. This another key ingredient of the calculation, and an accurate description of 
the infrared absorption can be only achieved using an accurate e-He interaction. We want 
to mention that DF theory has been succesfully applied to dynamical problems in bulk 
liquid helium and helium droplet o 22 ' 23 ' 24 using a generalization of the OP density functional 
called Orsay-Trento (OT) functional.— Having no ad-hoc parameters to be adjusted in the 
calculation to describe the infrared-absorption energies, we will show that the simplest form 
of a DF for liquid helium, namely a zero range one, is not able to reproduce the experimental 
data. As we are going to show, one must resort to more sophisticated forms, such as the 
mentioned OP and OT finite range functionals. This finding is an additional supporting 
evidence of the reliability of finite range functionals. 

Density functional theory has also proven to be the most successful approach in ad- 
dressing cavitation in liquid helium so far . 26 ' 27 ' 28 It incorporates in a self-consistent way the 
equation of state of bulk liquid and surface tension of the liquid-gas interface as a function 
of temperature. It allows for a flexible description of the electron bubble, incorporating sur- 
face thickness effects. Within DF theory one avoids the use of macroscopic concepts such as 
surface tension and pressure at a nanoscopic scale; however, it is a continuous, not an atomic 
description of the system. Used in conjunction with a Hartree-type e-He potential, we have 
recently shown that this approach quantitatively reproduces the existing experimental data 
on cavitation of electron bubbles in liquid helium below saturation pressure.— 9 Consequently, 
it is a tested framework to address other properties of electron bubbles. 

This paper is organized as follows. In Sec. II we briefly present the density functional 
(DF) plus Hartree electron-effective potential approach we have employed, as well as some 
numerical details. The structure of the electron bubble (radius and surface thickness) and 
infrared-absorption energies are discussed in Sec. Ill, and a Summary is presented in Sec. 
IV. Finally, we present in the Appendix a model calculation of the surface tension as a 
function of pressure that we have used to make some crude estimates of the electron bubble 
radius. 
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II. DENSITY-FUNCTIONAL APPROACH AND ELECTRON-HELIUM INTER- 
ACTION 



As in Ref. |29|, our starting point is the finite temperature zero-range DF of Ref. 30 that 
reproduces thermal properties of liquid 4 He such as the experimental isotherms and the 4 He 
liquid-gas coexistence line up to T = 4.5 K, and the T dependence of the surface tension of 
the liquid free surface. We have taken the Hartree-type e-He effective potential derived by 
Cheng et al.— (see also Ref. l32r ) as e-He interaction. This allows us to write the free energy 



of the system as a functional of the 4 He particle density p, the excess electron wave function 
tf, and T: 

F[p,vD,T] = J dff(p,T) + ^-/ df\VW)\ 2 + J dr\*(f)\ 2 V(p) , (2) 

where f(p,T) is the 4 He free energy density per unit volume written as 

f(p, T) = f vol (p, T) + + £(Vp) 2 . (3) 

In this expression, f vo i(p,T) consists of the free energy density of a Bose gas, plus phe- 
nomenological density and temperature dependent terms written as 

hnMT) = h(T)p>+ l -c(T)p^ (4) 

that takes into account the effective interaction of helium atoms in the bulk liquid. The 
parameters of these terms and those of the density gradient terms in Eq. © have been 
adjusted so as to reproduce physical quantities like the equation of state of the bulk liquid 
and the surface tension of the liquid free surface. 

The e-He interaction V(p) is written as a function of the local helium density 3 ^- 

, , Ti 2 kl 2-nh 2 „ /4vr\ 1/3 dn 

^=2^ + ^'°«- W (t) < 5) 

where a = 0.208 A 3 is the static polarizability of a 4 He atom, and ko is determined from 
the helium local Wigner-Seitz radius r s = (3/47TP) 1 / 3 by solving the trascendental equation 

t&n[k (r s - a c )} = k r s , (6) 

with a c and a a being the scattering lengths arising from a hard-core and from a polarization 
potential. We have taken 31 a a = —0.06 A, a c = 0.68 A. 



For given P and T values we have solved Euler-Lagrange equations which result from the 
variation of the constrained grand potential density u(p, ^, T) = u(p,^,T) — e\^\ 2 , where 

u{p, vl/, T) = f(p, T) + ^ |V* | 2 + M 2 V{p) - pp . (7) 

It yields 

ft 2 

' -Atf + V(p)V = , (9) 



2 m e 

where e is the lowest eigenvalue of the Schrodinger equation obeyed by the electron. These 
equations are solved assuming spherical symmetry, imposing for p that p'(0) = and p(r — > 
oo) = pb, where pb is the density of the bulk liquid, and that the electron is in the Is state. 
Fixing pb and T amounts to fix P and T, since the pressure can be obtained from the bulk 
equation of state P = —f vo i{Pb,T) + ppb, and the 4 He chemical p = df vo i(p,T)/dp\T is 
known in advance. We have used a multidimensional Newton- Raphson method^ for solving 
Eqs. (JHJ) and (|9*)l. after having discretized them using 13-point formulas for the r derivatives. 
A fine mesh of step Ar = 0.1 A has been employed, and the equations have been integrated 
up to Roo = 150 A to make sure that the asymptotic bulk liquid has been reached. After 
obtaining the equilibrium configuration (Is state), the spectrum of the electron bubble e n i 
is calculated from Eq. (JHJ) keeping frozen the helium density (Franck- Condon principle). 



III. RESULTS 
A. Liquid 4 He 

Figure H shows the Is — lp transition energies (eV) as a function of P (atm) for T = 1.25 
K. It can be see that the agreement between theory and experiment is good from P = 
to the solidification pressure, and consequently, the physical process seems well undestood. 
However, a minor discrepancy appears at high pressures, and especially in the description of 
the Is — 2p transition energies, as it can be seen in Fig. |21 This is not surprising, since one 
would expect that the 2p state is more sensitive to fine details of the bubble structure, in 
particular to its thickness, because the 2p wave function penetrates deeper into the liquid. 

We have improved the method of Sec. II to achieve a better agreement with experiment. 
The improvement is based on the observation that zero-range DF's as the one described 
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in Sec. II are fitted to reproduce the experimental surface tension of liquid helium at zero 
temperature , 30 ! 34 and then, the T-dependence of a,— and the thickness t of the free surface 
-defined as the difference between the distances at which the density equals 0.1p b and 0.9pb, 
where pt, is the bulk density at the given (P, T)- come out as predictions of the formalism. 
It turns out that at T = 0, zero-range DF's overestimate t by about 1 A. Indeed, recent 



measurement of the 4 He surface thickness yield values around 6 A, Refs. 



value predicted in Refs. 



30 



34 



35 



36, whereas the 



is about 7 A. In the case of 3 He, the experimental value is about 
7.5 A,— whereas the prediction using zero-range DF's is about 8.5-9 A . 34 i 38 We recall that 
the surface thickness of liquid helium is a quantity rather difficult to determine experimental 
and theoretically— and the dispersion of the values assigned to it is l ar g e . 34 i 35 i 36 i 37 i 38 i 39 i 40 i 41 
Only recently, the mentioned values (~ 6-6.5 A for 4 He and ~ 7.5 A for 3 He) have emerged 
as likely accurate determinations of the surface thickness of liquid helium free surface. 

In the nineties, a new class of DF's has appeared that retains some of the simplicity of 
the original zero-range DF's, incorporating finite-range effects that are absolutely necessary 
to address a wide class of physical phenomena, like elementary excitations in bulk liquid 
and large inhomogeneities caused by the presence of impurities in bulk liquid and droplets, 
and also by substrates. Two such DF's are the OP and OT functionals for 4 He already 
commented.—^ It is remarkable that they reproduce the experimental surface tension at 
T = without having imposed it in their construction, and also the surface thickness. In 
particular, OP yields a = 0.277 K A~ 2 , t = 5.8 A, and OT yields a = 0.272 K A~ 2 , t = 6 
A-2 - 2 ^ A recent diffusion Monte Carlo calculation yields a = 0.281(3) K A~ 2 , £=6.3(4) A.ii 
It is then quite natural to ask oneself if the remaining disagreement between theory and 
experiment shown in Fig. (Q) and especially in Fig. @ might arise from a worse description 
of the electron bubble surface in the case of zero-range DF's. Indeed, the error in the infrared 
energies has been estimated^ to be ~ 1%, smaller than the difference between experiment and 
zero-range DF calculations. To answer this question, we have repeated the calculations using 
the zero temperature OT functional -thermal effects on the excitation energies are expected 
to be small for T = 1.25 K, and this is what we have found, see below-, and also using a 
straightforward generalization at finite T of the OP functional, which is possible because of 
the similarities between OP and the zero-range DF of Sec II.— It is worth to know that a 
generalization of the OT density functional at finite temperatures is also available. 43 

Density profiles of electron bubble at different pressures and T = 1.25 K, obtained with 
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the OP and the zero-range functional, are shown in Fig. El and in Fig. 0] for the OT 
functional at T = K. The excess electron squared wave functions |\1/| 2 corresponding to 
the Is, lp, and 2p states are also represented. For the zero-range DF we also give the value 
of the surface thickness t, which goes from 6.1 A at P = atm to 4.3 A at P = 20 atm. 
In the case of OP, and especially of OT, it is difficult to define the thickness, and even the 
bubble radius, because of the density oscillations in the surface region (see also Ref. |22J). For 
reference, we indicate that the surface thickness of the electron bubble at T = 1.25 K and 
P = atm obtained with OP is 5 A, one A smaller that the value yielded by the zero-range 
DF. 

Figures and (J2J) also show the results obtained with OP and OT. It can be seen that 
in this case, the agreement between theory and experiment is excellent.— It is worthwhile 
to mention that Eloranta and Apkarian-2, have also obtained the absorption energies up to 
P ~ 12 atm using the OT functional and the pseudopotential of Jortner et al.^ 5 as e-He 
interaction. Their results for the Is— 2p transition energies are not in as good agreement with 
experiment as ours, likely because of the different e-He interaction we have used. Besides, 
the relevance of a correct description of the surface thickness to reproduce the experimental 
data has been overlooked in their work, which is mainly about the dynamics of electron 
bubble expansion in liquid 4 He. 

Figure (0) shows the evolution of the radius of the electron bubble with pressure deter- 
mined by several authors using different methods , 4 i 46 i 47 i 48 and the value we have obtained 
using the zero-range and the OP and OT functionals. For a given pressure, we have defined 
the radius of the bubble as the radial distance at which the helium density reaches the cor- 
responding bulk value divided by two. It can be seen that our results are compatible with 
the more recent experimental determinations, even those obtained with the zero-range DF, 
which overestimate the radius of the electron bubble by ~ 0.5 A. 

We present in Fig. [HI the variation of the Is — lp excitation energy with T at P = 2.9 
atm. Open circles are the experimental data of Grimes and Adams, Ref. U, and the dashed 
line is the result obtained by these authors using their Spherical-square-well (SS W) model 
with the T dependence of the surface tension as at saturation vapor pressure. 16 Our results, 
displayed as filled circles, have been obtained using the finite temperature OP functional. 
It can be seen that above T = 2 K we overestimate the excitation energies. In spite of the 
good agreement between experiment and the SS W model, we believe it is coincidental, and 
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that the failure of our method to describe the experiment above T = 2 K is the physical 
signature of the breaking down of the Franck-Condon principle at high T. Indeed, above 
T ~ 2 K, it is no longer valid to assume that the electron bubble is void of He atoms. It 
has been shown^^ that the bubble is progessively filled with helium vapor that must be 
'dragged' by the electron in the course of its dipole excitation. This dynamical effect has to 
be taken into account if one wants to have a quatitative agreement with experiment above 
T ~ 2 K. Although it is in principle possible to address this problem within time-dependen ; 



43, 



density functional theory 2 ^ using the finite temperature OP functional or that of Ref. 
this is not a trivial issue and it is beyond the scope of our work. 

Finally, we have also calculated the cross-section of the transition from the Is to the 
lp and 2p excited states. For photons of frequency w polarized along the z axis, the cross 
section is given by 

4vr 2 e 2 

a(w) = — — w\(np\z\ls)\ 2 5{E np - E ls - hw) (10) 

We collect in Table U the excitation energies (eV), oscillator strengths and total cross 
sections (eV cm 2 ) to the first two p states obtained with the OT functional at T = 0, P = 
atm. Our results are in full agreement with these of Ref.—. 



B. Liquid 3 He 

We have also studied the infrared-absorption spectrum for electron bubbles in liquid 3 He. 
The application of DF theory to describe electron bubble explosions in 3 He proceeds as 
shown in Sec. II. We have used the zero-range DF proposed in Ref. |38| to describe the 



inhomogeneous liquid, and also a finite-range DF, called from now on FR, obtained from 



20 for the 4 He case. The e-He 



the zero-range one following the procedure indicated in Ref. 
interaction is given again by Eq. (J3J with the parameters corresponding to 3 He, namely 
a = 0.206 A 3 , and same values for a a and a c . 

We represent in Fig. [7| the infrared-absorption energies for the Is — lp and Is — 2p 
transitions. It can be seen that the results are qualitatively similar to those found for 4 He, 
i.e., a fair insensitivity of the Is — lp transition energy to the detailed structure of the bubble 
surface, and an underestimation of the Is — 2p transition energy by the zero-range DF. The 
half density radius of the electron bubble K\n goes from 22.5 A at P = to 11.8 A at 
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P = 22.3 atm. 

The smaller excitation energies in the case of 3 He are due to the smaller surface tension 
for this isotope, a = 0.113 K A -2 , 49 which causes that electron bubble radii are larger for 
3 He than for 4 He (e.g., R eq = 23.5 A at P = instead of 18.9 A), yielding smaller excitation 
energies, as shown for instance by the simple model of Eq. In particular, at P = 

the model gives for the ratio of the Is — lp excitation energies the value (as/a^) 1 ^ 2 , in good 
agreement with the results obtained within the DF approach (see Tables). 

Finally, we collect in Table HP the excitation energies (eV), oscillator strengths and total 
cross sections (eV cm 2 ) to the first two p states obtained with the FR functional at T = 0, 
P = atm. It can be seen that the oscillator strengths and cross sections are sensibly the 
same for both isotopes. 

IV. SUMMARY 

We have demonstrated the suitability of the density functional approach to quantitatively 
address electron bubbles in liquid He. On the one hand, we have shown that the DF 
approach, in conjunction with a realistic electron-helium interaction, is able to reproduce 
without any further assumption, the low temperature infrared spectrum of electron bubbles 
experimentally determined by Grimes and Adams4^ On the other hand, the method yields 
results that agree with the experimental findings of Maris and co-workers on cavitation of 
electron bubbles below saturation pressure,— 1 ^ in a wide range of temperatures.— 

We have shown that the analysis of infrared-absorption transitions of electron bubbles 
constitute a stringent test for the theoretical models aiming a detailed description of the 
free surface of liquid helium, in particular of its surface thickness. This is at variance with 
electron bubble cavitation, which seems to be sensitive only to global properties of the 
surface, like its tension. Indeed, we have checked that the cavitation pressure for electron 
bubbles in liquid 4 He yielded by the zero-range DF is -2.07 bar, whereas OP yields -2.13 
bar, and OT yields -2.08 bar. These differences are far smaller than the experimental error 
bars, and are partly due to the small differences in the surface tensions predicted by these 
functionals. 

There are some related problems that can be studied as a natural extension of the work 
carried out until present within the DF frame. In particular, the effect on the critical 
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cavitation pressure of quantized vortices pinned to excess electrons.— The infrared absorption 
above T ~ 2 K is another open problem whose quantitative description requires to relax the 
Franck-Condon principle. Although they were not detected in the experiments,— ^ it would 
be also interesting to check if there are other resonances in the infrared spectrum arising 
from the state where the helium atoms are relaxed, even though this state itself is not 
realized. For instance, in the case of electrons in crystals, this phenomenon, pertaining to 
polaron physics, is the origin of a peak in the electron optical spectrum on top of the Franck- 
Condon spectrum, due to the resonance of this virtual relaxed excited state with the applied 
radiation field.— Both problems, presently unaffordable by any microscopic approach, can 
be addressed within DFT and TDDFT. As these are not trivial issues at all, they may be 
the next test grounds to assess the capabilities and limitations of the DF theory applied to 
liquid helium. 



APPENDIX 

In this Appendix we work out the zero temperature surface tension of bulk liquid 4 He as 
a function of P for the zero range DF. As we have previously stated, this is a model system, 
since bulk helium can only be in equilibrium with its vapor along the coexistence line, and 
thus for one single P value at the given T. 

We closely follow the method developed in Refs. Isflssl We take the planar interface 
perpendicular to the z axis, and consequently, the liquid density is only function of z. At 
given P, the equation of state determines the value of the bulk density p& and chemical 
potential \i. Starting from Eqs. (j3J) and (jlj) taken at T = 0, one writes the Euler-Lagrange 
equation 

d ^-p( P - + 2 p -)-2ip" = ^ (A.l) 
dp \ p p J 

where the prime denotes the z derivative. Eq. (jA.lj) admits solutions that go to p^ when 

re Euler-Laerange 



z — > — oo, and to zero when z — > oo. For t 

p(z) -actually z(p)- as indicated in Refs. 3011381 Density profiles corresponding to P = 0, 
5, 10, and 20 bar are shown in Fig. |H] as a function of z. Note that for each profile, the 
location of the z = point is arbitrary, as it stems from the structure of Eq. (jA.lj) . and 
that the surface region is steeper the higher the pressure is. The variation with pressure of 
the surface thickness is shown in Fig. El 
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To obtain the surface tension, we first define a reference system in which surface effects 
are not taken into account. It is a step-density system of bulk density pb filling the z < 
halfspace, whose free surface is arbitrarily located at z — 0. Placing -also arbitrarily- the 
overpressurized 4 He density profile in such a way that p(z = 0) = Pb/2, we obtain a(P) 
subtracting the grand potential of the overpressurized system, uo = f(p) — pp, from that of 
the reference system, u>b = fvoi(pb) — PPb, if z < 0, and anywhere else, and dividing the 
difference by the surface of the planar interface. 

Defining 

A/ = fvol(p) ~ fvol(pb) 

Ap = p- Pb (A.2) 

we have 



o-(P) = 2 f° (Af-pAp)dz 



oo 



+ 2 / [f vol (p) -p P + p^+ £p' 2 ] dz . (A.3) 



io p 
Using that Eq. ()A.1|) can be integrated to yield 38 



dp_^( A/ - P& P * ' " 
dz { £ + £ 

a close formula for the surface tension can be readily obtained: 



(A.4) 



a(P) = 2 f Pb (Af - pAp) 1 / 2 (P+A dp + 

Jpb/2 \P j 

Note that the explicit knowledge of p(z) is not needed to determine a. This is due to the 
particular form of Eq. (jA.lj) . 38 

Fig. shows the surface tension (K A~ 2 ) as as function of P (atm). It can be seen that 
in the displayed pressure range, <j(P) increases almost linearly. The variation we find from 
to 20 atm is similar to that found in Ref. Q. 
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FIG. 1: 4 He Is — lp transition energies (eV) as a function of P (atm) for T = 1.25 K. The open 
triangles are the observed points from the electron bubble photocurrent,— and the open circles 
correspond to direct infrared absorption measurements^ The dashed line represents the results 
obtained using the zero-range DF discussed in Sec. II. Filled circles connected with a dotted line 
are the results obtained using the OP finite-range DF at T = 1.25 K, and filled diamonds connected 
with a solid line are results obtained using the zero temperature OT finite-range DF. 

FIG. 2: 4 He Is — 2p transition energies (eV) as a function of P (atm) for T = 1.25 K. The 
open triangles are the observed points from the electron bubble photocurrenti 4 ' The dashed line 
represents the results obtained using the zero-range DF discussed in Sec. II. Filled circles connected 
with a dotted line are the results obtained using the OP finite-range DF at T = 1.25 K, and filled 
diamonds connected with a solid line are results obtained using the zero temperature OT finite- 
range DF. 

FIG. 4: Electron bubble density profiles in A~ 3 (right scale) and excess electron squared wave 
functions |^| 2 in A -3 (left scale), as a function of the radial distance r (A) for T = K and 
different pressures, obtained using the OT functional. The electron squared wave functions \^\ 2 
are represented as in Fig. |3J i.e., dotted line, Is state; dot-dashed line, lp state; thin solid line, 
2p state. The vertical thin lines indicate the equilibrium radius R eq yielded by the simple electron 
bubble model Eq. (JTJ using the values of cr(P) obtained in the Appendix. 



FIG. 3: Electron bubble density profiles in A~ 3 (right scale) and excess electron squared wave 
functions |^| 2 in A -3 (left scale), as a function of the radial distance r (A) for T = 1.25 K and 
different pressures. The dashed lines correspond to the zero-range DF, and the solid lines to the 
OP functional. The vertical thin lines indicate the equilibrium radius R eq yielded by the simple 
electron bubble model Eq. Q using the values of <r{P) obtained in the Appendix. The value of 
the surface thickness t is also given for the zero-range DF. The electron squared wave functions 
1^1 2 correspond to the OP calculation; dotted line, Is state; dot-dashed line, lp state; thin solid 
line, 2p state. 
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FIG. 5: Radius of the electron bubble as a function of pressure. Dashed line, zero-range DF 
results; filled circles(squares), OP(OT) results at T = 1.25(0) K; dotted line, Ref. U results. Open 



squares, data with error bars, and open triangles are from Refs. 



% respectively. 



FIG. 6: Variation of the Is — lp excitation energy with T at P = 2.9 atm. Open circles are the 



experimental data of Ref. |5J. Dashed line, SSW model results^ Filled circles, finite temperature 
OP functional results. The solid line has been drawn to guide the eye. 



FIG. 7: 3 He Is — lp and Is — 2p transition energies (eV) as a function of P (atm) for T = K. 



The dashed line is the result obtained using the zero-range DF of Ref. 



38. Filled dots connected 



with a solid line are the results obtained using the zero temperature finite-range DF. 

FIG. 8: Density profiles (A) for the overpressurized liquid 4 He model system described in the 
Appendix, as a function of r (A) for several values of P (bar). The results have been obtained 
using the zero range DF. 

FIG. 9: Surface tension (K A~ 2 , left scale) at zero temperature as a function of P (atm) for the 
overpressurized liquid 4 He model system described in the Appendix. Also shown is the surface 
thickness t (A, right scale). The results have been obtained using the zero range DF. 

TABLE I: Transition energies, oscillator strengths and total cross sections for the first two dipole 
excitations of electron bubbles in 4 He calculated with the OT functional at T = and P = atm. 



Transition 


Transition 


Oscillator 


Total cross 




energy (eV) 


strength 


section (eV cm 2 ) 


ls-lp 


0.105 


0.971 


1.07 xl0~ 16 


ls-2p 


0.488 


0.0250 


2.75 x!0~ 18 



TABLE II: Transition energies, oscillator strengths and total cross sections for the first two dipole 
excitations of electron bubbles in 3 He calculated with the FR functional at T = and P = atm. 



Transition 


Transition 


Oscillator 


Total cross 




energy (eV) 


strength 


section (eV cm 2 ) 


ls-lp 


0.071 


0.972 


1.07 xlO" 16 


ls-2p 


0.326 


0.0248 


2.73 x!0~ 18 
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